KrMTVIlT  or  CMUKL/NiIdmi  Bwiiii  (I  SUMiili 


DISCLAIMER  NOTIC 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


Reproduced  by  the 

CLEAR  INGHOUSE 

tor  fodoral  Scientific  &  Technical 
Information  Sprmgfiold  Va.  22151 


CONTENTS 


ABSTRACT . iii 

1.  INTRODUCTION . 1 

2.  PERVASIVE  DAMPING . 5 

3.  STABILITY  INEQUALITIES . 30 

REFERENCES . 41 

FIGURES 

1.  Isomorphic  System  -  Most  General  Conical  Case  ....  8 

2.  Stability  Diagram  -  Hyperbolic  Case 

(where  K  =  -0.  175  and  =  2.  34) . 36 

3.  Orientation  of  the  Principal  Axes  of  Inertia  of  the 

Composite  Satellite  System  -  Hyperbolic  Case . 37 

4.  Stability  Diagram  -  Hyperbolic  Case 

(where  K  =  -0.  165  and  *  =  2.  34) . 38 

5.  Stability  Diagram  -  Hyperbolic  Case 

(where  K  =  -0.175,  and  *>  =  2.  44) . 39 

TABLES 

1.  Conditions  for  Pervasive  Damping . 24 


v 


Air  Force  Report  No. 
SAMSO-TR-69-430 


Aerospace  Report  No. 
TR-0066(5306)-8 


PERVASIVE  DAMPING  AND  STABILITY  WITH  APPLICATION  TO 
A  GRAVITATIONALLY  STABILIZED  SATELLITE  AUGMENTED. 
BY  A  CONSTANT  SPEED  ROTOR 


Prepared  by 

R.  K.  Williamson 
Electronics  Division 


69  JUL  15 


Engineering  Science  Operations 
THE  AEROSPACE  CORPORATION 


Prepared  by 

SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
AIR  FORCE  SYSTEMS  COMMAND 
LOS  ANGELES  AIR  FORCE  STATION 
Los  Angeles,  California 


This  document  has  been  approved  for  public 
release  and  sale;  its  distribution  is  unlimited 


FOREWORD 


This  report  is  published  by  The  Aerospace  Corporation,  El  Segundo, 
California,  under  Air  Force  Con'  :act  No.  F04701  -69-C-0066. 

This  report,  which  documents  research  carried  out  from  September  1968 
to  December  1968,  was  submitted  01  69  Nov  07  to  SAMSO  (SMTO)  for  review  and 
approval. 


Approved 


A.  J.  ^chiewe,  Director 

Control  and  Sensor  Systems  Subdivision 

Electronics  Division 

Engineering  Science  Operations 


Publication  of  this  report  does  not  constitute  Air  Force  approval  of  the  report's 
findings  or  conclusions.  It  is  published  only  for  the  exchange  and  stimulation 
of  ideas. 


si  3  _ 

'jack  S.  Friedman,  1st  Lt.,  USAF 
‘Project  Officer 


ii 


ABSTRACT 


The  efficacy  of  a  particular  satellite  design  cannot  be  guaranteed 
without  first  examining  the  problem  of  pervasiveness  of  damping.  Two 
recent  contributions  in  the  field  of  motion  stability  analysis  to  the  general 
problem  of  identifying  whether  a  system  is  pervasively  damped  are 
discussed.  If  the  damping  is  indeed  pervasive,  then  a  corollary  to  the 
Thomson- Tait-Chetaev  stability  theorem  can  be  used  to  establish  condi¬ 
tions  for  which  the  system  is  asymptotically  stable.  These  results  are 
illustrated  in  terms  of  a  gravitationally  stabilized  vehicle  that  has  been 
augmented  by  means  of  constant  speed  rotors. 


1. 


INTRODUCTION 


This  report  considers  the  stability  of  a  gravitationally 
stabilized  satellite  that  has  been  augmented  by  means  of  a  constant 
speed  rotor. 

The  essential  features  of  the  strong  form  of  the  Thomson-Tait- 
Chetaev  theorem  are  reviewed,  and  its  relevance  to  problems  in  atti¬ 
tude  dynamics  is  discussed.  The  notion  of  pervasive  damping  is 
treated,  along  with  a  detailed  exposition  of  two  recent  papers  in  the 
literature  that  specifically  treat  the  problem  of  identifying  systems 
that  are  pervasively  damped. 

The  derivation  of  the  linearized  differential  equations  of 
small  angle  motion  about  equilibrium  can  be  found  in  Ref.  1. 
Generically,  they  are  of  the  form 


Mx  +  Gx  +  Dx  +  Kx  =  0 


(1) 


where  (as  was  defined  in  Ref.  1) 

x  is  a  state  vector  of  displacementt- from  equilibrium 
M  is  a  symmetric  positive  definite  inertia  matrix 
G  is  a  skew- symmetric  matrix  arising  from  gyroscopic  forces 
D  is  a  symmetric  matrix  arising  from  damping  in  the  system 
K  is  a  symmetric  matrix 

The  strong  form  of  the  Thomson-Tait-Chetaev  theorem  can 
immediately  be  deduced  from  a  standard  application  involving  a 
Lyapunov  function.  The  function  V  is  said  to  be  a  Lyapunov  function 
if  the  following  conditions  are  fulfilled: 

a.  V(x,t)  is  defined  in  ft  for  all  t  >  0  where  ft  is  a  set, 
defined  as  ft:  i  x  ||<  A  and  t  t  0,  over  which  the  existence 
and  unicity  theorem  holds  for  the  system  of  differential 
equations. 

b.  V(0.t)  =  0  for  t  >  0 
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c.  V(x,t)  dominates  a  certain  u(x)  for  all  x  in  ft  and  all 
t  >  0 

d.  V(x,t)<0infl 

The  system  under  consideration  is 


x  =  f (x,  x,t) 

The  interpretation  of  x  in  the  original  definition  of  a  Lyapunov  function, 
however,  can  be  broadened  to  include  this  special  system.  Consider 
the  Lyapunov  function 

V=  1/2xTMx+  1/2xTKx  (2) 

The  time  rate  of  change  of  V  is 


^j—  =  xTMx  +  xTKx 


If  Eq.  (1)  is  premultiplied  by  x  ,  then 


x^Mx  +  x^Dx  +  x^Gx  +  x^Kx  =  0 


Since  G  was  assumed  to  be  skew-symmetric,  the  following  identity  can 
be  obtained: 


xTMx  +  xTKx  =  -xTDx 


The  time  rate  of  change  of  the  Lyapunov  function  is  given  by 


dV  =  _*TD* 
dt  x  x 


The  properties  of  the  matrices  that  are  necessary  to  insure  that  the 
equilibrium  solution  is  asymptotically  stable  can  immediately  be 
deduced  from  an  application  of  a  theorem  originated  by  Lyapunov. 
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Under  the  conditions  of  the  definition  of  a  Lyapunov  function,  it 
can  be  further  established  that  there  exists  a  continuous,  non¬ 
decreasing  scalar  function  p(x)  such  that  j3( 0)  * **  0  and,  for  all  t, 

V(x,  t)  <  p(  i'x  ||) 

and  V(x,  t)  is  strictly  decreasing  for  all  solutions  0  except  the  null 
solution,  then  the  equilibrium  solution  is  asymtotically  stable  and  any 
other  solutions  in  0  converge  to  the  equilibrium  solution.  Upon  appli¬ 
cation  of  the  theorem  wc  find  that:  (1)  if  M,  K,  and  D  are  positive, 
definite,  V  is  positive  definite,  V  is  negative  definite,  and  asymptotic 
stability  is  assured:  and,  (2)  due  to  an  instability  theorem  originally 
due  to  Lyapunov,  if  M  or  K  is  negative  definite  or  sign  variable  and 
D  is  positive  definite,  V  is  not  pngitiv-  definite,  V  is  negative  definite, 
and  the  equilibrium  solution  is  unstable  Thus,  for  a  positive  definite 
damping  matrix,  the  stability  or  instaoility  of  the  equilibrium  solution 
can  be  inferred  from  the  inertia  and  stiffness  matrix  alone.  Since  the 
inertia  matrix  is  guaranteed  to  be  positive  definite,  asymptotic  stability 
or  instability  can  be  determined  solely  from  the  stiffness  matrix  alone. 
Collectively,  these  results  constitute  the  "strong -form"  of  the 
Thomson-Tait-Chetaev  theorem. 

In  the  absence  of  damping,  it  is  possible  to  have  a  stable 
equilibrium  solution,  but  examination  of  the  eigenvalues  of  K  alone  do 
not  lead  to  all  possibilities  for  stability  l  2  J  .  In  general,  other  possi¬ 
bilities  exist  in  which  the  G  matrix  plays  a  role.  Systems  that  rely 
upon  the  gyroscopic  matrix  for  stabilization  will  be  unstable  upon  the 

>5*  v 

introduction  of  the  appropriate  damping  into  the  system. 


*  Numbers  in  brackets  refer  to  references. 

**  Appropriate  damping,  at  least  at  this  juncture,  means  damping  such 
that  a  positive  definite  damping  matrix  exists. 
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The  assumption  of  a  positive  definite  damping  matrix  is  too 
restrictive.  In  the  design  of  attitude  control  systems,  one  strives  not 
for  a  positive  definite  damping  matrix,  but  rather  for  damping  that 
affects  the  entire  system,  so  that  any  deviation  from  the  nominal 
motion  elicits  energy  dissipation.  Systems  for  which  energy  dissipa¬ 
tion  exists  for  any  motion  of  the  system  are  said  to  be  pervasively 
damped.  The  criteria  of  pervasive  damping  for  systems  whose  pri¬ 
mary  mechanism  for  the  dissipation  of  energy  is  via  viscous  damping 
are 


Dx  <  0 


for 


♦  (t.XQ.to) 


and 


T 

-x  Dx  =  0  for  x  =  0 

where  Mt.Xg.tg)  represents  a  possible  motion  of  the  system.  This 
definition  does  not  necessarily  entail  the  imposition  of  a  positive 
definite  damping  matrix.  At  this  juncture,  we  are  in  a  position  to  state 
a  corollory  to  the  "strong-form"  of  the  Thomson-Tait-Chetaev 
theorem.  If  we  can  construct  linearized  differential  equations  of  small 
angle  motion  about  equilibrium  in  the  form  suggested  in  Eq.  (1).  with 
a  positive  definite  inertia  matrix,  a  gyroscopic  matrix,  a  symmetric 
matrix  K.  and  a  damping  matrix  derivable  from  a  Raleigh  dissipation 
function  that  is  not  necessarily  positive  definite,  and  if  further  the 
damping  is  pervasive,  then  the  equilibrium  solution  is  asymptotically 
stable  if  the  matrix  K  is  positive  definite. 

These  results  serve  to  motivate  some  basic  ideas  encapsuled  in 
a  theorem  that  is  of  paramount  utility  in  motion  stability  analyses.  The 
theorem  was  originally  advanced  by  R.  Pringle  [  3  j  and  subsequently 
refined  by  P.  W.  Likins  f  4  J  and  by  P.  W.  Likins  and  D.  L.  Mingori  [  5  ]. 

Theorem 

A  holonomic  mechanical  system  with  the  Lagrangian  not 
explicitly  dependent  upon  time  and  with  damping  both  solely  generalized 
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and  pervasive  in  the  noncyclic  generalized  coordinates  is:  (1) 
asymptotically  stable  in  these  coordinates  if  the  dynamic  potential 
energy  is  positive  definite,  and  (2)  unstable  in  these  coordinates  if  the 
dynamic  potential  energy  P  can  take  on  negative  values  arbitrarily 
close  to  the  origin  in  the  coordinate  space. 

In  essence,  the  Hamiltonian  is  advanced  as  an  appropriate 
Lyapunov  function.  The  Hamiltonian,  which  is  represented  in  terms  of 
the  potential  energy  and  homogeneous  forms  determined  from  the 
kinetic  energy,  is  given  as 


H  =  T2  +  v  -  T0 

=  T2  +  P  (4) 

Since  the  homogeneous  form  T2  can  be  guaranteed  to  be  positive 

definite,  the  positive  definiteness  of  the  Hamiltonian  is  established 

solely  by  examination  of  the  dynamic  potential  energy  P. 

The  link  between  the  two  separate  approaches  relies  upon  the 

observation  that  the  Lyapunov  function  in  Eq.  (2)  is  a  quadratic 

approximation  of  the  Hamiltonian  [  4  J.  In  addition,  the  inner  product 
T 

x  Kx  is  identified  as  a  first  approximation  to  the  dynamic  potential 
energy.  The  Hamiltonian  {or  its  approximation)  is  positive  definite  if 
the  dynamic  potential  energy  (or  its  approximation)  is  positive  definite. 

2.  PERVASIVE  DAMPING 

In  most  practical  applications,  the  damping  matrix  of  Eq.  (1) 
is  not  positive  definite.  We  have  seen  that  this  provides  no  real 
impediment  to  the  determination  of  whether  an  equilibrium  solution  is 
asymptotically  stable  or  unstable,  providing  the  damping  is  pervasive. 
Determination  of  whether  a  system  is  pervasively  damped  or  not  was 
at  best  uncertain  prior  to  the  recent  publication  of  two  articles  that 
specifically  deal  with  this  problem.  R.  E.  Roberson  establishes  a 
constructive  method,  using  the  concepts  of  graph  theory,  for  deter¬ 
mining  from  the  structure  of  the  dynamical  equations  alone  if  a  linear 
system  is  pervasively  damped  [6  ].  G.  M.  Connell  subsequently 
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advanced  an  alternate  method  [7  ].  Essentially,  one  determines  a 
manifold  of  solutions  for  which  there  exists  no  damping,  and  proceeds 
to  show  that  the  manifold  of  solutions  do  not,  in  fact,  satisfy  the  dif¬ 
ferential  equations.  Thus,  any  possible  system  motions  elicit  energy 
dissipation.  These  two  alternative  methods  of  determining  whether  a 
dynamic  system  is  pervasively  damped  will  be  elaborated  upon  in  more 
detail. 

R.  E.  Roberson  starts  with  the  observation  that  the  system  of 
equations  represented  in  Eq.  (1)  is  isomorphic  to  a  set  of  n  point 
masses  free  to  move  in  one  dimension  and  interconnected  by  lumped, 
linear,  passive,  mechanical  elements,  such  as  springs,  viscous 
dampers,  and  gyrators.  The  linear  differential  equation  of  motion  for 
the  point  mass  system  in  matrix  form  is 


ml 

0 

0 

0 

0 
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This  can  be  conveniently  accomplished  by  an  orthogonal  transformation 
of  the  original  equations  of  motion.  Defining  the  orthonormal  trans¬ 
formation  such  that  the  M  matrix  is  diagonalized,  we  have 

Mz+Dz+Gz+  Kz  =  0 


where 


z  =  Px 

The  identification  can  be  made  complete  by  selecting  the  masses  of  the 
isomorphic  point  mass  systems  equal  to  the  diagonal  elements  of  M'. 
The  matrix  sum  D'  +  G'  is  then  rewritten  as  a  sum  of  a  symmetric  and 
a  skew- symmetric  matrix.  Next,  the  point  mass  isomorphic  system 
is  represented  in  a  diagram.  The  most  general  representation  for  the 
conical  case  appears  in  Fig.  1.  It  is  assumed  that  the  configuration 
parameters  are  consistent  with  the  general  representation  that  appears 
in  Fig.  1.  An  investigation  in  terms  of  a  diagram  is  then  made  to 
determine  whether  the  system  is  completely  connected  or  whether  it 
separates  into  several  disconnected  pieces.  For  the  latter  alternative, 
damping  in  one  part  of  the  system  will  not  be  experienced  in  another 
portion  of  the  system.  To  facilitate  this  investigation,  a  connective 
matrix  U  is  defined  as  follows.  The  connective  matrix  U  has  as  many 
rows  as  there  are  point  masses  comprising  the  system,  and  it  has  as 
many  columns  as  there  are  mechanical  element  connections  between 
the  point  masses.  If  a  mechanical  element  j  is  connected  to  a  point 
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Figure  1.  Isomorphic  System  - 

Most  General  Conical  Case 

mass  i,  then  a  1  is  entered  in  that  matrix  element  (i,  j)  corresponding 
to  the  ith  point  mass  and  the  jth  connecting  element.  In  this  way,  a 
connection  matrix  can  be  evaluated.  For  the  case  examined  here,  the 
matrix  U  ia 


11  10  0  0 


U  = 


1 

0 


0 


0 

1 

0 


0  110 
0  10  1 
10  11 


If  the  rows  and  columns  of  U  can  be  rearranged  in  any  way  such  that  it 
taxes  the  form 


U 


* 

1 


0 


0 


U 


* 


0 


0 


0 


0 
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that  is,  of  a  diagonal  partitioned  matrix  with  zero  off-diagonal  blocks, 
then  the  system  contains  as  many  separate  parts  as  there  are  sub¬ 
matrices  on  the  diagonal  of  the  partitioned  matrix.  Roberson  develops 
an  algorithm  that  is  beneficial  for  assessing  the  number  of  partitions 
in  the  event  that  the  system  is  extremely  complex.  The  existence  of 
separated  subsystems  assures  that  the  damping  is  non-pervasive. 
Roberson  observed  that  a  sufficient  condition  for  instability  is  that  the 
columns  of  the  damping  and  stiffness  matrices  are  linearly  dependent. 
If  a  motion  exists  such  that  there  is  no  dissipation  of  energy,  the 
coordinates  must  all  remain  equal.  Moreover,  any  coordinate  cor¬ 
responding  to  a  particle  connected  by  a  damper  to  inertial  space  must 
remain  zero.  Denote  a  matrix  V  consisting  of  those  columns  of  U 
associated  with  a  dissipative  element.  The  matrix  V  is  subjected  to 
the  same  analysis  that  the  matrix  U  was  subjected  to.  For  the  case 
examined,  wc  find 


Let  Ta  denote  those  point  mass  elements  associated  with  the  ath  sub¬ 
graph  of  the  matrix  V.  The  set  T  is  comprised  of  exactly  one  element 
from  each  set  TQ.  For  the  case  considered  here,  we  define 


T:M2 

Tj  :  Mg,  M^i 


Define  also  a  set  Z  comprised  of  the  union  of  the  zero  rows  of  matrix 
V  and  set  T.  The  set  Zq  denotes  those  indices  corresponding  to  a 
particle  connected  by  a  damper  to  inertial  space.  Thus,  Z  and  ZQ  are 
comprised  of  the  elements 
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Z: 

Zq!  M2,M3,  m4 

The  complete  motion  is  described  by  a  new  set  of  coordinates  .  .  . 
u>r  smaller  in  number  than  the  Z.  and  specifically  equal  to  the  number 
of  elements  in  set  Z  diminished  by  the  nu  iber  of  elements  in  Zq. 

The  transformation 


Z  =  Qu 


where  Q  is  an  n  X  r  matrix  constricted  algorithmically  as  follows.  If 
i<Z0’  the  ith  row  of  Q  is  a  row  of  zeros.  If  i/(TaUZQ),  enter  a  digit  1 
in  the  ith  row,  ath  column  of  Q  and  fill  the  remainder  of  the  row  with 
zeros.  Otherwise,  if  u(T0UZq),  enter  a  digit  1  in  the  ith  row  at  the 
first  column  that  does  not  already  contain  a  digit  1  in  some  other  row. 
This  specific  process  has  the  effect  of  assigning  label  numbers  to  the 
u>.  which  describe  the  motion  of  the  a  part  of  the  matrix  V.  The 
remaining  point  mass  indices  on  u.  comprise  those  masses  which  are 
neither  a  part  of  V  nor  are  fixed  in  position  with  respect  to  inertial 
space.  The  resulting  matrix  is 


1 


Q  = 


0 

0 


One  can  observe  from  Fig.  1  that  in  order  for  there  to  be  no  dissi¬ 
pation  of  energy,  point  masses  2,  3,  and  4  cannot  move  relative  to 
each  other,  and  further,  there  can  be  no  motion  of  the  subpart  com¬ 
prised  of  the  point  masses  2,  3,  and  4.  When  transformed,  the 
equations  of  motion  become 

mQu  +  gQiL  +  (k  +  K)  Qu  -  0 
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aMMS*******  -***■ 


Substituting  the  trial  solution 


<*>  =  exp(Xt)wQ 

we  obtain 

WUo=  0 

where 

W  =  [x2m  +  Xg  +  (k  +  K)]Q 

A  necessary  and  sufficient  condition  that  a  non-trivial  solution  exists 
is  that  the  rank  of  W  is  not  maximum.  The  rank  can  be  checked 
algorithmically  as  follows.  Choose  any  r  rows  of  W  and  determine 
the  X  values  that  make  its  determinant  vanish.  Having  X  equal  to  one 
of  these  r  values  is  a  necessary  condition  for  rank  less  than  r.  There 
remain  n!/r!(n-r)!  -1  possible  determinants,  each  to  be  evaluated 
r  times.  If  one  encounters  a  non-zero  determinant  at  any  point  in  the 
evaluation  sequence,  it  is  guaranteed  that  no  non-trivial  solution 
exists.  If  the  rank  of  W  is  not  maximum,  then  there  exists  a  system 
motion  for  which  there  is  no  dissipation  of  energy,  and  then  the  system 
is  not  pervasively  damped.  For  the  case  considered  here, 


If  gj,  g2»  and  g3  =  0,  X  is  identically  zero.  Letting,  X  =  0  in  the  first 
equation,  w  will  be  of  maximum  rank  if  K  4  0.  Therefore,  one  might 
conclude  the  system  will  be  pervasively  damped  for  those  systems  for 


which  the  isomorphic  representation  is  as  described  in  Fig.  1.  Or, 
2 

solving  for  X  from  the  first  equation  (K  i  0)  and  substituting  the 
2 

values  of  X  into  the  equations 


x2g2  *  o 
x2g2  4  0 
x2g3  4  0 

we  can  verify  that  W  is  of  maximum  rank  if  gj,  g^,  or  g^  is  not 
equal  to  zero. 

An  alternate  approach  to  the  determination  of  whether  a  system 
is  pervasively  damped  is  given  presently  [  7  ].  This  method  was  sug¬ 
gested  by  G.  M.  Connell.  The  power  dissipation  function,  or  the  V  in 
Eq.  (3),  is  given  as 


F(x)  =  -xTDx 

where  the  damping  matrix  for  the  conical  and  hyperbolic  case  is  of  the 
form 


[o  0  0  0 


D  - 


0 

0 


Lp 


0 

0 

0 


0 

0 

0 


0 

0 

D 


We  wish  to  determine  those  solutions  for  which  there  is  no  dissipation 
of  energy.  A  necessary  and  sufficient  condition  for  there  to  be  no 
energy  dissipation  is  that  the  control  body  be  fixed  relative  to  the  main 
body.  This  condition  is  represented  by  specifying  that  the  state  vari¬ 
able  component  x^  be  some  arbitrary  constant.  There  exist  eight 
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possible  conjectured  system  motions  for  which  there  is  no  dissipation 
of  energy,  and  they  are 


where  c.  represents  a  state  variable  that  is  considered  to  be  some 
arbitrary  constant  other  than  zero,  and  denotes  a  state  variable  that 
is  not  considered  constant. 

If  this  manifold  of  solutions  for  which  there  exists  no  damping 
does  not  satisfy  the  differential  equations,  then  the  damping  is  per- 
vasive.  Consider  the  solution  vector 


Substituting  the  solution  back  into  the  differential  equations  of  motion, 
we  have 


A .  v* 

A  t  + 


.* 


(6) 


where 


i 


* 


x2 


and  the  matrices  A  ,  G ' ,  and  K  are  [n  x  (n  -  m)] . 

% 

The  essence  of  the  following  amounts  to  the  determination  of 
those  conditions  for  which  the  only  possible  solution  of  equations  is  the 
null  solution.  Thus,  the  manifold  of  solutions  for  which  there  exists 
no  damping  is  not  satisfied  by  the  differential  equations.  In  pursuit  of 
these  objectives,  consider  the  homogeneous  part  of  Eq.  (6). 

Assuming  the  trial  solution 


C  =  exp{\t)$J 

and  substituting  into  the  homogeneous  equation  we  find 

W  C  -  0 


where 


W  =  X2A'  +  \G'  +  K 
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.  ittj. 


In  total,  there  exist 


n! 

(n  -  m)!m! 

different  set*  of  X's  obtained  by  selecting  m  rows  of  W  and  finding 
those  values  of  X.  for  which  the  determinant  vanishes.  A  necessary  and 
sufficient  condition  for  a  non-trivial  solution  to  the  homogeneous  equa¬ 
tion  to  exist  is  that  the  rank  of  W  be  less  than  maximum.  If  the  rank 
of  any  one  of  the  possible  W's  is  maximum,  then  the  only  solution  that 
exists  is  the  trivial  one.  On  the  other  hand,  if  after  the  evaluation  of 
the  different  W's,  it  is  determined  that  no  non-zero  determinants 
exist,  the  rank  would  not  be  maximum  and  thus  one  is  guaranteed  that 
a  non-trivial  solution,  other  than  the  equilibrium  solution,  exists. 

Thus,  damping  would  not  be  pervasive  in  this  instance.  The  particular 
solution  of  Eq.  (6)  is  found  from 

K'C  =  C' 


since  C '  is  constant. 

Consider  the  augmented  [nX  (n  ■  m  +  l)j  matrix  determined  by 
combining  the  (n  -  m)  columns  of  the  matrix  K'  with  the  column 
matrix  C'.  A  necessary  and  sufficient  condition  for  a  non-trivial 
solution  to  exist  is  that  the  rank  of  the  matrix  K'  is  equal  to  the  rank 
of  the  augmented  matrix  (K't  C').  If  it  is  determined  that  the  rank  of 
both  matrices  is  equal,  then  a  non-trivial  solution  exists. 

Having  established  the  methodology,  we  will  proceed  to  apply  it 
to  the  hyperbolic  case.  The  equations  of  motion  for  the  specialized 
hyperbolic  case  were  generically  of  the  form 
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Consider  the  possible  solution 


from  the  set  of  eight  possible  solutions  for  which  there  exists  no 
damping.  First,  let  us  examine  the  particular  solution 

Kl*  =  C' 

where 


For  the  hyperbolic  case  and  the  possible  solution  under 
examination 


"a 

b 

0' 

V 

b 

c 

0 

K2{ 

K  = 

0 

0 

d 

.  C'  = 

0 

e 

f 

0 

*2* 

where  K 2  is  an  irbitrary  constant.  If  for  this  particular  solution  one 
can  assume  that  the  K'  and  the  augmented  matrices  are  of  maximum 
rank,  then  the  trivial  solution  is  the  only  one  possible.  Since  ele¬ 
mental  operations  on  a  matrix,  such  as  multiplication  of  one  column 
by  a  constant,  do  not  change  the  rank  of  the  matrix,  one  can  see  that 
the  augmented  matrix  is  of  maximum  rank  if  the  stiffness  matrix  K  is 
of  the  maximum  rank.  This  is  assured  if  the  stability  inequalities  are 
satisfied.  The  K'  matrix  is  also  guaranteed  to  be  of  maximum  rank  if 
three  of  the  stability  inequalities  are  satisfied.  The  conclusion  that 
the  only  solution  is  the  trivial  one  depends  upon  the  examination  of  the 
homogeneous  solution. 

The  W  matrix  for  the  hyperbolic  case  and  the  particular  solution 
examined  is 


0 

0  " 

"0 

0 

-*1 

a 

b 

0‘ 

1 

0 

m2 

0 

0 

0 

"g2 

b 

c 

0 

3 

(1 

0 

0 

m3 

+  \ 

gl 

g2 

0 

+ 

0 

0 

d 

1 

_0 

m4 

0  _ 

0 

to 

0 

g3 

e 

f 

0 

In  total,  there  exist  four  X.  sets  obtained  from  the  solutions  of  the 
four  possible  determinants,  A  necessary  and  sufficient  condition  for 
a  non-trivial  solution  to  exist  is  that  the  rank  of  W,  for  all  solutions 
of  X  obtained  from  the  four  different  \  sets,  be  less  than  maximum. 
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In  the  event  that  the  rank  of  W  is  maximum  for  any  solution  X,  the  only 
solution  that  exists  is  the  trivial  one.  The  four  equations  obtained  by 
equating  these  four  possible  determinants  to  zero  are 


^(mjm^mj)  +  K^arn^m^  +  cm^m^  +  dm^m^  +  gZm2  + 

2  2  2  2 
+  X  (ac  +  arr^d  +  cdnrij  -  2bgjg2  +  gjC  -  b  +  ag2) 


+  (acd  -  b  a)  =  0 


(7) 


4  2 

X  (g3m1m2  +  g2m1m4)  +  X  (g3am2  +  g-jcrnj  -  gjtt^b 
+  gjem2  +  g2am4  +  g^nij)  +  (ag;}c  -  bg2e  -  gjfb 

+  gjec  -  bZg3  +  g2af)  =  0  (8) 


4 

X  (em2m3  -  g2g1m4  -  gjg3m2  -  bm^) 

2  2 
+  *■  (bg2g3  +  cm3e  +  cdm2  -  gjg2f  +  g2e 

"  cgjg3  "  bfm3  -  bdm4)  +  (ecd  -  bfd)  =  0  (9) 

\6(-m1m4m3)  +  X^g^nij  -  gjim4  -  m3am4  -  fm^j 

2  2 

-  dmjm4)  +  X  (ag2g3  +  bem3  -  gjf  +  g1g2e 

-  bgjg3  "  afm3  -  adm4  -  fdmj)  +  (bde  -  daf)  =0  (10) 

Solutions  for  X  identically  equal  to  zero  have  been  specifically 
excluded.  The  particular  degenerate  case  is  considered  later.  For  a 
non-trivial  solution  to  exist,  all  the  determinants  would  have  to  be 
identically  zero  for  all  solutions  of  X.  Thus,  Eq.  (7)  would  have  to 
be  satisfied  by  the  solution  to  Eqs.  (8),  (9),  and  (10),  and  similarly, 
for  the  remaining  equations. 


1  8 


|  46  * 

s 


ml 


It  can  be  observed  immediately  that  Eqs.  (8)  and  (9)  will  be 
satisfied  by  the  roots  to  Eqs.  (7)  and  (10),  only  if  the  latter  two 
equations  have  repeated  roots.  Further,  it  will  be  observed  that  the 
roots  of  each  equation  must  span  the  space  of  roots  corresponding  to 
the  other  equations  in  order  for  all  determinants  to  vanish  identically, 
and  that  W  will  be  less  than  full  rank.  The  examination  of  seven  addi¬ 
tional  possible  system  motions  for  which  there  exists  no  damping 
remains.  One  additional  case  will  be  elaborated  in  the  text,  then  a 
table  that  summarizes  calculations  associated  with  the  remaining 
possible  motions  will  be  presented. 

Consider  the  possible  solution 

xl" 
c2 
X3 
C4 

from  the  remaining  set  of  possible  motions  for  which  there  exists  no 
damping.  The  matrices  K'  and  C'  become 


a  O' 

bkj  +  ek2 

b  0 

ckl  +  fk2 

c'  = 

0  d 

0 

e  0 

/kl  +  gk2. 

where  Kj  and  K2  are  arbitrary  constants.  The  matrix  K'  is  guaranteed 
to  be  of  full  rank  if  the  stability  inequalities  are  satisfied.  Specifically, 


i 


a 


v 


The  augmented  matrix  is 


0  bkj  +  ek2 

0  ckj  +  fk2 

d  0 

0  fkj  +  gk^ 

The  rank  of  the  augmented  matrix  can  be  guaranteed  to  be 
maximum  if  one  of  the  four  possible  determinants  obtained  by  alter¬ 
nately  deleting  a  row  is  non-zero.  One  of  the  determinants  was  iden¬ 
tically  zero  since  the  matrix  contained  a  column  of  zeros.  The  three 
remaining  determinants  are 


kj(af  -  eb)  +  k^(ag  -  e^)  =  0 

HI) 

kjd(b^  -  ac)  +  k2d(eb  -  af)  =  0 

(12) 

kjd(bf  -  ec)  +  k2d(bg  -  ef)  =  0 

(13) 

Consider  a  solution  to  Eqs.  (11)  and  (12)  for  the  constants  kj 
and  k^,. 


"af  -  eb 

2  i 

ag  -  e 

r-i .  n 

d(b^  -  ac) 

d(eb  -  af)_ 

5V  " 
tN> 

II 

'  Q- 

The  solution  for  the  two  constants  can  be  guaranteed  to  be  other  than  a 
trivial  one  if  the  determinant  is  identically  zero,  that  is 


d(af  -  eb)^  +  d(ag  -  e^)(b^  -  ac)  =  0 
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If  the  stability  inequalities  are  satisfied,  then  it  can  be  verified  that 
the  determinant  is  not  equal  to  zero,  and  the  only  solution  is  the  trivial 
one.  If  kj  and  kg  are  equal  to  zero,  then  the  homogeneous  solution  is 
the  only  one  to  be  considered.  Two  additional  cases  remain.  Consider 
a  solution  to  Eqs.  (12)  and  (13)  for  the  constants  kj  and  kg: 


d(b^  -  ac) 

d(eb  -  af) 

V 

o' 

d(bf  -  ec) 

d(bg  -  ef^ 

1 

.k2_ 

0 

Again,  the  solution  for  the  two  constants  can  be  guaranteed  to  be  other 
than  a  trivial  one  if  the  determinant  is  identically  zero.  Since  d  is 
guaranteed  to  be  greater  than  zero, 

(b^  -  ac)(bg  -  ef)  -  (eb  -  af)(bf  -  ec)  =  0 


or 


b(acg  +  2bfe  -  e^c  -  b^g  -  af^)  =  0 

Since  the  term  in  parentheses  is  guaranteed  to  be  greater  than  zero, 
the  determinant  is  identically  zero  if  b  is  equal  to  zero.  So,  if  b  is 
equal  to  zero,  a  solution  for  kj  and  k g  exists  other  than  the  trivial  one. 
In  fact,  if  b  is  equal  to  zero,  kj  and  kg  are  related  by 


If  this  relation  is  substituted  into  Eq.  (11),  we  find 

t2  ,  2?  n 
-  af  +  cag  -  ce  =0 

If  the  relation  above  is  identically  zero,  then  we  are  guaranteed  that  the 
rank  of  the  augmented  matrix  will  be  less  than  maximum. 
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Similarly,  solving  for  kj  and  k2  from  Eqs.  (11)  and  (13)  we 
find  that  if  e  is  equal  to  zero,  and  if  further 

(b2  -  ac)g  +  af2  =  0 

then  the  rank  of  the  augmented  matrix  is  not  maximum.  Now,  if  both 
e  and  b  are  identically  equal  to  zero,  Eqs.  (11)  and  (12)  become, 
since  a  is  greater  than  zero, 


kjf  +  k2g  =  0 
kjC  +  k2f  =  0 

and  Eq.  (13)  vanishes  identically.  A  solution  other  than  the  trivial 
one  exists  if  the  determinant  is  equal  to  zero.  Thus,  for  the  case  when 
both  b  and  e  are  equal  to  zero,  if  further 

f2  -  eg  =  0 

then  the  rank  of  the  augmented  matrix  is  not  maximum.  Now,  the 
homogeneous  solution  will  be  considered.  The  matrix  W  for  the 
solution  considered  is 


“ml 
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-«l" 

a 

o' 

1 

0 

0 

0 

-«2 

b 

0 

3- 

ii 

* 

0 

m3 

+  X 

*1 

0 

+ 

0 

d 
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0  _ 

0 

g3. 

e 

0 

22 


Consider  the  six  possible  determinants  which  have  been  equated  to 
zero,  as  follows: 


-  ggitijX2  +  gjb  "  g2a  s  0  <14) 

m^m^X4  +  (am3  +  nijd  +  g2)X2  +  ad  s  0  (15) 

g3mjX2  +  g3a  +  gje  *  0  (16) 

(bm3  +  gjg2)^2  +  bd  *  0  (17) 

(*>g3  +  g2®)k=  0  0®) 

(gjg3  -  em3)X2  -  ed  =  0  (19) 


If  the  solution  to  any  one  of  these  six  possible  equations  is  not  satisfied 
by  any  of  the  remaining  solutions,  then  the  rank  of  the  W  matrix  is 
maximum.  Thus,  the  only  system  motion  would  be  the  trivial  one. 

For  example,  if 


*>g3  +  g2®  4  0 


from  Eq.  (18),  then  the  determinant  is  not  equal  to  zero  for  any 
solution  of  X  other  than  a  trivial  one.  The  trivial  solution  is  ruled  out 
since  it  implies  that  one  of  the  variables  is  constant,  and  those  cases 
are  specifically  considered  separately.  The  calculations  associated 
with  the  remaining  six  possible  motions  as  well  as  the  calculations 
already  presented  in  the  text  are  enumerated  in  Table  1. 

We  are  left  with  eight  specialized  system  motions  for  which  there 
is  no  damping.  Additional  conditions  can  be  delineated  such  that  damp¬ 
ing  exists  only  for  the  trivial  solution.  Damping  exists  for  solutions 
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Table  1.  Conditions  for  Pervasive  Damping 


Table  1.  Conditions  for  Pervasive  Damping 


Case  Homogeneous  Solution  Particular  Solution 
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Table  1.  Conditions  for  Pervasive  Damping  (Cont) 


if  b,  f  =  0  and  e 


other  than  the  trivial  one,  and  hence  is  not  pervasive,  if  simultaneously 
e,  f  and  g  are  zero  or  if  simultaneously  bg  -  ef  =  0,  bf  -  ce  =  0  and 


A  limiting  example  of  the  hyperbolic  case  is  where  the  constant 
speed  rotor  is  oriented  along  the  pitch  axis.  The  hinge  axis  is  also 
oriented  along  the  pitch  axis.  Librations  in  roll  or  yaw  would  not  be 
accompanied  by  energy  dissipation,  and  so  this  particular  configuration 
would  not  be  pervasively  damped.  This  will  be  confirmed  presently  in 
terms  of  the  preceding  equations.  The  equations  of  motion  for  this 
special  case  reduce  to 
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Consider  the  possible  system  motion 

x'  =  [xt  c2  x3  c4]T 


Consider  the  six  possible  determinants  which  have  been  equated  to 
zero,  Eqs.  (14)  through  (19).  If  the  solution  to  only  one  of  these 
six  possible  equations  is  not  satisfied  by  any  of  the  remaining  solutions, 
then  the  rank  of  the  W  matrix  is  maximum,  and  a  non-trivial  solution 
would  exist.  For  the  special  case  examined,  all  but  one  equation 
remains,  since  the  parameters  g^,  g^,  b,  and  e  are  identically  zero. 
Thus,  the  rank  of  W  is  not  maximum,  and  we  can  conclude  that  the 
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configuration  for  the  special  case  of  pitch  augmentation  with  a  hinge 
axis  along  the  pitch  direction  would  not  be  pervasively  damped. 

To  this  juncture,  we  have  considered  two  different  methods  for 
determining  whether  a  system  is  pervasively  damped.  Definitive 
judgments  about  specific  configurations,  i.  e.  ,  configurations  for 
which  the  parameter  values  are  known  explicitly,  can  be  made.  More 
general  statements  about  the  pervasiveness  of  damping  for  classes  of 
configuration  parameters  is  more  difficult,  although  the  relationships 
that  must  be  satisfied  in  order  for  there  to  be  no  pervasive  damping 
are  enumerated. 

3.  STABILITY  INEQUALITIES 

In  developing  the  stability  inequalities,  we  find  that  the  existence 
of  a  positive  definite  stiffness  matrix  is  sufficient  to  insure  that  the 
system  is  asymptotically  stable,  if  the  damping  is  pervasive.  The 
stiffness  matrix  for  the  hyperbolic  case  is  represented  generically  as 


K  = 


a 

b 

0 

e 


b  0  e“ 
c  0  f 
0  d  0 
f  0  g 


The  K  matrix  can  be  guaranteed  to  be  positive  definite  if  the  following 
inequalities  are  satisfied: 


a  >  0  (20) 

ac  -  b^  >  0  (2 1) 

d  >  0  (22) 

aeg  +  2bfe  -  e^c  -  b^g  -  f^a  >  0  (23) 
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For  the  conical  case,  the  stiffness  matrix  is  generically  represented 
by 


fa  0  0  0 


K  = 


0 

0 


b 

e 


IP  t 


e  f 

c  g 

g  d 


The  stability  inequalities  for  the  conical  case  are 


a  >  0 

(24) 

b  >  0 

(25) 

be  -  e2  >  0 

(26) 

bed  +  2efg  -  f2c  -  e2d  -  bg2  >  0 

(27) 

The  stability  inequalities  associated  with  the  conical  case  will  not  be 
developed  further,  since  the  performance  of  this  particular  configura¬ 
tion  is  markedly  inferior  to  the  performance  experienced  with  the 
hyperbolic  case.  This  comparison  is  made  more  precise  in  Ref.  8. 

For  the  hyperbolic  case,  the  relationship  between  the  symbols 
employed  in  the  stiffness  matrix  and  the  system  parameters  is  as 
follows : 


a 


(‘I  '  «3 
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b  = 
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J0  (i*  - 1.)  _ 

=  "  j^-Sv?sg  -  — j- — !  (3  +  s2e) 

/  T  \ 

=  — -  (c2e  -  s2e)  +-=£  (cec^-  s*>se) 

1  1 


e  =  sees. 


f  =  -  ~  (3  +  s2e) 
1 


g  =  -  y-  (3  +  S2e)  +  £ 

A1  A1 


Consider  the  following  definitions: 


k 


A1!- 


T 


I 


1 

2 


1 

1 


k 


These  definitions  are  not  identical  to  the  definitions  for  inert: 
parameters  that  are  commonly  used  in  dynamic  analyses. 


(28) 


32 


Since  the  main  body  of  the  spacecraft  is  symmetric,  and  the  control 
body  is  a  rod  with  a  hinge  axis  along  the  body  pitch  axis,  the  following 
relationships  hold: 


Also 


The  symbols  employed  in  the  stiffness  matrix  reduce  to 


(29) 


a  =  k4c?ce  -  kj  (c2e  +  3) 

b  =  kjSGCG  -  k^C^SG 

c  =  -  kjSpSG  -  (kj  +  k3)(s26  +  3) 

d  =  k3(c26  -  S2e)  +  k4(C?C6  -  S<pS0) 
e  =  k3S6C6 

f  =  -  k3(s2e  +  3) 
g.-k3(s2e  +  3)  +  ^. 
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As  defined  in  Ref.  1,  the  equation  that  must  be  satisfied  for  equilibrium  is 


sin  (6  +  V)  =  K  sin  26  30) 

So,  if  K  and  <p  are  given,  there  exist  several  6's  that  satisfy  the 
equation  above, but,  in  general,  only  one  value  of  6  gives  results  that 
are  stable.  This  will  be  considered  in  more  detail  later; 

Considering  the  multidimensionality  of  the  parameter  space,  it 
is  difficult  to  give  a  complete  representation  of  the  stability  in  terms 
of  all  the  configuration  parameters.  The  size  of  the  parameter  space 
could  be  reduced  in  part  by  limiting  the  values  that  the  parameters  are 
allowed  to  take  based  upon  the  performance  analysis  results. 

If  the  parameters  <p  and  K  were  given,  there  would  exist  at  least  one 
value  of  6  from  Eq.  (30).  The  proper  selection  of  6  depends  upon 
satisfaction  of  the  inequality  in  Eq,  (22).  Since  K  is  specified  in 
advance,  there  exists  a  relationship  between  the  parameters  and  k^, 
where  k^  is  a  measure  of  the  inertia  properties  of  the  control  body 
relative  to  the  main  body  roll  axis  and  k^  is  a  measure  of  the  rotor 
"strength"  relative  to  the  inertia  of  the  main  body  about  the  roll  body 
axis.  The  inequality  in  Eq.  (22)  becomes 

(A  -  sze)  -  (cjc^fse)  >  0 


or 


CO.  28  -  co"jg+  6>  >  0  (31) 

So,  given  a  value  of  K  and  9,  there  exists  a  value  of  6  that  satisfies 
Eq.  (30)  and  the  inequality  in  Eq.  (31).  Since  the  parameters  k^ 
and  k^  are  linearly  related  as  in  Eq.  (29),  there  exist  three  inde¬ 
pendent  quantities,  kj,  k^,  and  S/I,  which  serve  to  define  the  config¬ 
uration  uniquely.  By  selecting  the  parameter  K,  we  are  constraining 
rotor  "strength"  to  be  linearly  related  to  the  parameter  k^.  By 
changing  the  value  of  the  parameter  k^  we  can  change  the  rotor  strength 
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relative  to  a  fixed  value  of  the  parameter  k^,  but  we  also  have  a 
different  equilibrium  orientation  as  implied  by  Eq.  (30).  The  values 
of  K  and  are  based  upon  the  performance  results  in  Ref.  8.  In 
Fig.  2,  a  value  of  k  equal  to  -0.  175  and  *  equal  to  2.  34  radians 
was  selected.  The  value  of  G  was  3.  775  radians.  The  orientation  of 
the  body  is  illustrated  in  Fig.  3.  For  a  given  point  in  k^,  k ,  space, 
values  of  the  dimensionless  parameter  S/I  and  values  that  are  greater 
identify  the  stable  region. 

In  Fig.  4  the  parameter  value  of  k  was  changed  to  -0.  165. 
implying  a  greater  rotor  "strength"  for  a  comparable  value  of  the 
parameter  k^.  But,  the  equilibrium  position  is  changed  from  a  value 
of  6  equal  to  3.  775  radians  to  a  value  of  6  equal  to  3.  784  radians  as 
well.  For  a  given  value  of  k^,  and  the  increased  rotor  "strength,  " 
the  instability  boundary  associated  with  the  dimensionless  parameter 
S/I  shows  an  almost  imperceptible  change.  In  Fig.  5,  the  value  of 
was  changed  to  2.  44  radians,  and  the  value  of  K  was  equal  to  -0.  175. 
The  equilibrium  orientation  would  necessarily  change  as  well.  Again, 
the  instability  boundary  associated  with  the  dimensionless  parameter 
S/I  shows  a  very  slight  change  for  a  given  parameter  k^. 
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Stability  Diagram  -  Hyperbolic  Case  (where  K  =  -0.175  and  -  2.34) 
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